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Cumulant ratios in fully developed turbulence 
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In the context of random multiplicative cascade processes, we derive analytical solutions for one- and two- 
point cumulants with restored translational invariance. On taking ratios of cumulants in Ine, geometrical effects 
due to spatial averaging cancel out. These ratios can successfully distinguish between splitting functions while 
multifractal scaling exponents and multiplier distributions cannot. 



1. Introduction 

The Navier-Stokes equation governing fluid 
flow is deterministic. Nevertheless, the statisti- 
cal description of fully developed turbulence has 
a long tradition ||^. Random multiplicative cas- 
cade models form a particularly simple and ro- 
bust class of such statistical models, reproduc- 
ing important observed features such as mul- 
tiplier distributions and their correlations 
and related Kramers-Moyal coefficients reflect- 
ing Markovian properties The models have 
worked almost too well in the sense that different 
cascade-generating probability density functions 
(pdf's or "splitting functions") p{qL,qR} for the 
multiplicative weights have been equally success- 
ful in reproducing these observables. More so- 
phisticated ways to distinguish between them are 
clearly desirable. 

While some experiments have concentrated on 
measuring statistics in the energy dissipation den- 
sity £, we recently found a complete analytical 
solution working in In e rather than e itself . 
We here and in Ref. j|] show that cumulants in 
Ine are analytically calculable even when trans- 
lational invariance is restored in order to emulate 
the spatial homogeneity of experimental turbu- 
lence statistics. Both one- and two-point cumu- 



lants turn out to be powerful tools which for third 
and fourth order differ not only in magnitude but 
even in sign for splitting functions which are indis- 
tinguishable in terms of other observables. Unlike 
multifractal scaling exponents, for example, such 
cumulants can be expected to distinguish between 
different models for sufficiently large experimen- 
tal samples. 



2. Analytical solution for random multi- 
plicative cascades 

Energy flux densities e are generated in the 
simplest multiplicative cascade models as fol- 
lows. In successive steps j — 1,...,J, the in- 
tegral scale L is divided into equal intervals of 
length Ij — = and dyadic addresses 

K = (ki ■ ■ ■ kj) with ki = or 1. At each step 
j, the energy flux density Sk-i-.-kj generates fluxes 
multiplicatively in the two subintervals via 



^ki---kjkj^i — Qki-'-kjkj^i ^ki---kj j 



(1) 
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where the random variables ql = qkf kjO and 
QB. = qki---kji for the left and right subintervals 
are drawn from a given cascade-generating pdf 
p{qLt<1r), independently of other branches and 
generations of the dyadic tree. When after J cas- 
cade steps the smallest scale ij = Ij = L/2'-' is 
reached, the local amplitudes of the flux density 
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field 



n iki-k, 



(2) 



at positions < t{K)={\ + X^/^i kj2-^~^) < 2'^ in 
units of T] are interpreted as the energy dissipation 
amplitudes which are to be compared to experi- 
mental time series converted to one-dimensional 
spatial series by Taylor's frozen flow hypothesis. 

We have shown previously that, since the 
product of multiplicative weights (||) becomes ad- 
ditive on taking the logarithm, 
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Figure 1. The hierarchical structure of the cas- 
cade is reflected in a corresponding structure in 
the source parameters A. 



\nekik2---kj = ^\nqk,...k, 



(3) 



the multivariate cumulant generating function for 
In e has the analytical solution 



lnZ(A, 



0---0, 



,Ai...i) 



(4) 



= In/expl ^ Afci...fc, lnefej...fc^ 



_,...,kj=0 

J 1 

j = l fci,...,fcj_i=0 

where the branching cumulant generating func- 
tion Q has arguments 



A 



E 



A 



ki ■■■kj 



(5) 



,fcj=0 



(see Figure 1) and is defined by the Mellin trans- 
form of the splitting function. 



g(AL,Afl) =ln 



j dqLdqRp{qL,qR)ql'' q^" 



,(6) 



Because of the simplicity of (^, Q can often be 
found analytically. 

A host of analytical predictions for statis- 
tics in Ine follow, starting with any and all n- 
multivariate cumulants obtained directly from 
In Z through 

C(ki,K2,- ••,«;„) = ((IneKi) • • • (lneK„))c (7) 

d" In Z 



for arbitrary dyadic bin addresses ki,...,k„. 
These multivariate cumulants in Ine are easily 
calculated since, due to the additivity of In Z in 
(^), they are simple sums of same-lineage cu- 
mulants Cn and splitting cumulants Cr,s in In q (see 
eqs. (|l§ and (|2|) below). 



((In 9) 
9"Q 



{{InqLYilnqn) 



(8) 



(9) 



where without loss of generality we have assumed 
Q{Xl, Xr) to be symmetric in its arguments. 
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Figure 2. Translational averaging using the mov- 
ing window technique. 



A=0 



3. Restoring translational invariance 

Before the above theoretical cumulants can be 
compared to experimentaUy measured ones, the 
issue of translational invariance must be dealt 
with. Clearly, the generating function ^ and 
its cumulants (|^) are not translationally invari- 
ant, in conflict with the homogeneous statistics 
characterising experimental results. Spatially ho- 
mogeneous statistics can, however, be emulated 
by creating a theoretical time series consisting of 
a chain of m adjacent independent cascade fields 
with 2"^ finest-scale bins each Q. In analogy to 
the experimental situation, an observation win- 
dow of width 2'^?7 is successively moved over this 
series in bin-sized steps, t = 1, . . . , (m— 1)2'^, suc- 
cessively "seeing" parts of adjacent cascade con- 
figurations: see Figure 2. 

A translationally invariant one-point moment 
density would thus be constructed as 

^ (m-l)2-' 

Pn = 7 -TTTJ i^^^^tT , (10) 

m^oo TO — 1 2'' — ' 
^ ' t=l 

which should be comparable to the experimen- 
tal time series. Operationally, this can be imple- 
mented by keeping only one cascade while aver- 
aging over many different cascade configurations, 
i.e. 

M 1 

^» = mE((i-^*)")-mE^"W' (11) 
t=i t=i 

with M = 2'^ = L/rj and ( ) denoting configura- 
tion averaging. Likewise, a translationally invari- 
ant two-point density with constant distance rj d 
(with d — 1,2,3,...) between the two bins would 
be simulated by two adjacent cascades, 

I ^' 

with pr^s{t,t+d) = {{lnetY{lnet+dy)- As shown 
in Figure 3, bin t+d at some stage in the sum- 
mation exceeds M = 2^^ and hence would refer to 
the right-hand cascade while t would refer to the 
left-hand one. 

Given that the model provides solutions in 
terms of cumulants, it is tempting to apply this 
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Figure 3. Translational averaging for two-point 
statistics. Indices t and t+d run over all bins with 
the "euclidean distance" d kept fixed. 



averaging prescription directly to cumulants also, 
i.e. to define 

CrAd) = j-iY.^rAt,t+d), (13) 

i = l 

using for Cr,s{t, t+d) theoretical expressions ob- 
tained from (|7|). However, experimental cumu- 
lants are derived from measured moments rather 
than the other way round (for example C2 = 
P2 - pI, C3 = p3 - 3p2 + 2p?), so that averaging 
over moments rather than cumulants is manda- 
tory for theory also. The proper procedure is 
hence to convert theoretical cumulants (|^) to mo- 
ments, average these over t, and then convert 
these back to cumulants for experimental com- 
parison. 

For the one-point case, this convoluted route 
becomes simple: the n-th order one-point cumu- 
lant Cn{t) = C{ki= . . . =Kn=t), given by 

C„(t) = Jc„, (14) 

is independent of position t so that translational 
averaging is trivial. The only remaining compli- 
cation is that J is not an experimental observable, 
and this is easily addressed by looking at cumu- 
lant ratios C„/C„' for which the J-dependence 
cancels. Ratios of translationally averaged one- 
point cumulants, 

£n ^ ^ ((lng)")e 

Cn' C„' ((Ingf)" )c 

being independent of J, should hence be directly 
comparable to experiment. 
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To demonstrate the quality of these cumulant 
ratios, we consider three model distributions, all 
with factorised splitting function 



(16) 



namely a binomial distribution (often also termed 
the "a model"). 



p{q) = 



OL2 



ai + 012 
ai + a2 



5(g-(l-ai)) 
6iq-{l + a2)) , 



(17) 



with parameters ai = 0.3 and a2 = 0.65, a log- 
normal distribution 



piq) 



1 



r2/K aq 



exp 



1 ^ 



(18) 



with parameter a = 

r(/3i + /?2) 



piq) 



r(/3i)r(/32) 



0.42, and a beta distribution 

g/3i~i(8-g)&"i(i9) 



with parameters /3i = 4.88 = /32/7 and q S (0, 8). 
The beta model is particularly appealing because 
it parallels the experimental situation where en- 
ergy conservation in three dimensions results in 
a non-energy-conserving one-dimensional projec- 
tion. Parameter values quoted are the result of re- 
quiring (q) = 1 and best fits needed to reproduce 
observed multiplier statistics which hence 
cannot distinguish between these three models. 

As shown in Figure 4, all three splitting func- 
tions also have almost identical multifractal scal- 
ing exponents r(n) = ln((7")/ln2. Since (g) = 1 
by construction, r(l) is zero for all three dis- 
tributions. For 71 = 2 we get t(2) = 0.26 for 
the first two distributions and 0.23 for the beta 
distribution, indistinguishable within the uncer- 
tainty of the experimental intermittency expo- 
nent At = 0.25 ± 0.05 [|ll[. We secondly note 
that even for n > 3 the r(n)'s for the binomial 
and the beta distributions remain indistinguish- 
able. Thirdly, all three distributions have a pos- 
itive skewness {{q — (g))^)/((g — (9))^)^^^ when 
measured in q and reproduce the observed multi- 
plier statistics, including their correlations. This 
has been shown for the binomial and log-normal 
in Ref. |Q. Numerical analysis of the beta dis- 
tribution yields similar results. Fits to observed 




4 5 6 7 

order n 



Figure 4. Multifractal scaling exponents T(n) 
for the binomial (circles, eq. (|l7|)), lognormal 
(squares, eq. (|l8|)) and beta (diamonds, eq. (ITT 
distributions. 



scaling exponents have not been performed be- 
cause it is not straightforward to compare theo- 
retical and experimental scaling exponents due to 
the finiteness of the inertial range |l2j . 

The above observables thus fail manifestly to 
distinguish between the different model distri- 
butions. By contrast, Figure 5 demonstrates 
that, while ci and C2 are almost identical for the 
three models, higher-order cumulants and cumu- 
lant ratios (|l^) are very different. For example, 
C3 = 0.05, 0.00 and —0.05 for the distributions 
(p7|), ( p^ and ( p^ respectively, so that the theo- 
retical cumulant ratios 



£3 
C2 



((Ing)^), 



((lng)2)^ 



(20) 



0.31 
0.00 
-0.25 



(binomial) 

(lognormal) 

(beta) 



lead to results that differ even in sign. If the 
present model assumptions are adequate, this 
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Figure 5. Same-lin eag e cumulants c„ for the bi- 
nomial (circles, eq. (|l7|)), lognormal (squares, eq. 
and beta (diamonds, eq. (|l9|)) distributions. 



sign difference should be seen in the experimental 
ratio 



)bs 



c. 



)bs 



((ln£)3) -3((ln£)2)(ln£)-F2(ln£) 
((lne)2) - (ln£:)2 



.(21) 



In fourth order, the same-lineage cumulant C4 has 
a different sign for the binomial and beta distri- 
bution, so that the ratios C4/C2 = C4/C2 of two- 
point cumulants come with a different sign, too. 
For the log-normal distribution, these ratios are 
of course again zero. 

4. Two-point cumulants and geometry 

Having demonstrated the advantages of mea- 
suring ratios of one-point cumulants, we now 
consider equivalent two-point ratios. The the- 
oretical two-point cumulant for two bins is 
found in terms of their mutual ultrametric dis- 
tance D > 0. As illustrated 0| in Figure 6, 
two bins Ki — {ki ■ ■ ■ kjkjj^i ■ ■ ■ kj) and K2 — 
{ki ■ ■ ■ kjk'j^^ ■ ■ ■ k'j), with fcj+i 7^ ^j+ii are sepa- 




Figure 6. Some examples for the ultrametric 
distance D in a cascade with J = 3. 



rated by an ultrametric distance D = J — j . 

Cr,s {D)^[J~ D)Cr+s +Cr,s- (22) 

Again we must restore translational invariance us- 
ing ( |l^ ) for theoretical densities p^.s rather than 
( |T^ ) for cumulants. In doing so, we must note 
that, since the density Pr,s{t, t+d) factorises when 
bins t and t + d belong to independent cascades, 
i.e. whenever M — d+l<t< M, the averaged 
moment splits up, 



1 



M-d 



f=l 

M 

E 

t=M-d+l 



1 



Pr{t)ps{t + d) 



(23) 



Analytic expressions for pr^s{t,t + d) are again 
readily derived by inserting the cumulants (|T^ ) 
and ( p2|) into the usual relations between n- 
variate moments and cumulants fictl and thence 
into (H). 

Since the cumulants c„ and c^.s are indepen- 
dent of t, this procedure clearly involves summa- 
tion of (J — D) over t as in (|2^) to create "geo- 
metrical coefficients" of c„ and Cr,s of the type 



M-d 



Gn(J,d) = ^ X {J 



D{t,t + d)y 



(24) 



with n = 0, 1, 2, . . ., where the dependence of the 
ultrametric distance D on the bin positions is 
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made explicit. These coefEcients are best eval- 
uated by changing the index of of summation, 



Gn{J,d) 



J 

E 

D=l 



p{D\J,d) (J-D) 



(25) 



with p{D\J,d) the (normalised) histogram func- 
tion counting the number of times D appears 
while t runs over its allowed values. Empirically, 
we find 

r {1<D<A) 
piD\J,d) = \ l-{d/2^) {D = A) (26) 

[ d/2° {A<D<J) , 

where A = \\0g2 d] is the ceiling of log2 d. Inser- 
tion of (|2^) into (^5|) leads to analytical expres- 
sions for the geometrical coefficients 

Go{J,d) = (1-2-Jd), (27) 

Gi (J, d) = iJ-A)~ 2d{2-^ - 2"^) , (28) 

G2{J,d) = {J - Af ~ AdiJ - A)2-^ 

+ 6d{2-^ - 2--') , (29) 

which in turn yield analytical results for the aver- 
aged two-point densities J)^ ^{d) of (p3|). Spatially 
homogeneous two-point cumulants are then con- 
structed via the inversion formulae llOtl 



(30) 
(31) 



Ci.iid) = Pi,M)-pI, 
C2,iid) = P2a(c^) - 2piPi i(d) 

- P2P1 + 2-pl , 
Cz.i{d) = P3,iW - SpiPs^lc^) - P3P1 

- 3p2Pi,i('^) + 6p?Pi,i(rf) 

+ &-p{p\ - Q-Pl , (32) 
C2,2(d) = P2,2(rf)-4piP2,i(£^)-2(Pi,i(d))' 

- P2 + 8p?Pi,i(d) 

+ 4p2Pl - Gpt ■ (33) 
With Eqs. (||)-(|9|), we arrive for s=l at 

CrAd)^Gi{J,d)Cr + l+Go{J,d)Cr,l. (34) 

This turns out to be equivalent to direct trans- 
lational averaging of cumulants (|l3|). For 57^ 1, 
however, such direct averaging is wrong and the 
full conversion from cumulant to moment to av- 
eraged moment and back to averaged cumulant is 



unavoidable. For r=s=2 we get, for example, 
C2,2(rf) = G2{J,d)cl + Gi{J,d) {a + 4.0201.1) 

+ Go(J,d)(c2,2 + 2c2_i) 

- 2[Gi(J,d)c2 + Go(J,d)ci,i]' , (35) 

where the additional terms are a consequence of 
the third term in the expression ( |3^ ) for G2,2('i)- 
Averaged two-point cumulants Cr.s of higher or- 
der r, s>2 exhibit similar structures. Translation- 
ally averaged n-point cumulants Cmi,...,m„ can be 
calculated by the procedure sketched above. 

Figure 7 shows explicit examples for Go and 
Gi for a cascade of length J — 5. As expected, 
Go reflects the trivial dependence of the splitting 
cumulant on the sum limits 1 < ^ < M—d, apart 
from the point d~0 for which no splitting cumu- 
lant enters at all. More interesting is the coeffi- 
cient for the same-lineage cumulant, Gi: it con- 
sists of a series of straight-line segments, chang- 
ing slope whenever d mod 2 = and ending at 
zero for d > 2'^/2. Since the form of Gi changes 
whenever d is a power of 2, approximate expo- 
nential behaviour of Cr,i as a function of d is to 
be expected; this is shown in Figure 8. The ex- 
ponential form would, however, be destroyed by 
any sizeable contribution of Go entering via the 
splitting cumulant, especially at larger d. 

The form of Gi can be further understood 
by considering an alternative formulation for 
p{D\J,d) in terms oi k = U0S2 "^J ' 



piD\J,d) 



(1 - d2-'')5j^D 



(36) 



,7-1 

E 



e{J-j-k) (l-2^~"'d) 

x(^,/-D-j - Sj-D-j+l) , 



(with 6n — 6n.o the Kronecker delta and 0{n) = 
whenever n < and 1 otherwise) since, with 
Gi^j:'L=iiJ-D)p{D\J,d), we find 



j-i 

Gi(j,d) = ^e(2'^- 



d){l~2'-^d), (37) 



which is a sum of straight-line contributions kick- 
ing in whenever d becomes smaller than 2'^~^ ,j = 
1,2,... This means that whenever d becomes 
smaller than some dyadic fraction of M, the two 
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Figure 7. Geometrical coefficients Gi and Go, Figure 8. Geometrical coefficients Gi and Gq, 

on linear scale for J— 5. on logarithmic scale for J=5 and J— 10, showing 

the approximately exponential behaviour of Gi . 



bins t and t+d can fall within the same {J—j)- 
scale subcascade so that Gi picks up new contri- 
butions from this scale. 

Translationally invariant cumulants Cr,s are 
constructed from these factors according to eqs. 
(|^-(|3|). Figure 9 shows by example Gi.i for 
the binomial (a model) and the corresponding 
energy-conserving p-model with pdf 

p{qL,qR) = [^S(qL - I - a) + ^5{qL - I + a)] 
X 5{qL +qR-2) , 

setting (for purposes of comparison) ai — a2 — 
a = 0.4. The a-model has c^.s = and hence 
contains no contribution from Gq but only from 
Gi and C2 = i [ln((l + a)/(l - a))]^ = 0.1795, 
while for the p-model both ci i = — C2 and C2 
contribute in (^4|). The resulting p-model Gi^i 
has the same peak at d=0 as the a-model but 
exhibits the familiar anticorrelation (negative cu- 
mulant) at larger d . Whether and for what d 
the Gi_i is negative depends, however, on the sum 
of same-side and splitting cumulant contributions 
rather than on the splitting cumulant alone. 



We further note that all models whose split- 
ting function factorises have zero translationally 
invariant two-point cumulants for d > 2''/2. 
Roughly, this can be translated into the state- 
ment that deviations of two-point cumulants from 
zero for "long" distances d (compared to an ad- 
mittedly fluctuating cascade size which we have 
modelled as a constant 2"^) would signal the non- 
factorisation of the splitting function and vice 
versa. 

Returning to cumulant ratios, we focus on cu- 
mulants Cr.ijd). If the splitting function fac- 
torises as in (uw, then the splitting cumulant Cr.s 
is zero and the two-point cumulant for d > I, 
becomes directly proportional to the geometrical 
coefficient Gi(J, d), 

CrAd)=Cr+lGi{J,d). (38) 
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Figure 9. Ci^i two-point cumulants for the bi- 
nomial and p- models with J = 5, showing con- 
tributions from same-lineage and splitting parts 
of the respective splitting functions. Solid lines 
represent the (scaled) geometrical coefficients. 



Taking ratios of two-point cumulants of different 
orders, 



CrAd) _ cr+i _ ((lng)'-+i) 



+1\ 



(39) 



the geometrical coefHcient drops out, so that 
these ratios become independent of d. This is 
an important observation as it grants access to 
properties of the pdf (cascade generator) even af- 
ter spatial homogeneity has been restored. Also, 
the d-independence of these ratios constitutes a 
severe test of the model assumptions entering the 
cascade models. Furthermore, the factorisation 
assumption can be tested since one- and two- 



point ratios are equal if ( [iq ) holds, 

Crs{d) Cr+l Cr+l 



Cr',l{d) Cr' + l 



(40) 



In this case, C2,i/Ci,i would assume the same 
numerical values as those in Eq. (po|), with simi- 
lar powers to discriminate between models. Two- 
point cumulant ratios would hence also predict 
clearly different results, in contrast with scaling 
exponents and multiplier distributions. 

We also note that the connection 1^ between 
the multifractal scaling exponents r(n=A) and 
the cumulant branching generating function (0), 



'(A) = 



In(g^) _ Q(A,0) 



In 2 



In 2 



(41) 



implies that the same-lineage cumulants 
which are to be extracted from ratios (^5|) or i 
are related to the T(n) by 



9"Q(A,0) 



5A" 



= ln2 



5"t(A) 



A=0 



(42) 



A=0 



i.e. the cumulant in Inq is related to the n-th 
derivative of the scaling exponent t(A), taken at 
A = 0. In principle, this not only allows for an 
unambiguous, albeit indirect extraction of scal- 
ing exponents, but also of the more fundamental 
splitting function. 

We end on an interesting sideline regarding the 
detection of scaling in the bin size i. Convention- 
ally, this is done by plotting ln(e") against \ni 
in the expectation of seeing a straight line. The 
same scaling of (e") is just as easily detected by 
pointing out that the one-point cumulant in Ine 
is given at every scale j by C^^^ = jc„ and, since 



In 2 



(43) 



in other words, scaling in (e") is manifest in a 
logarithmic dependence on the length scale £j 
of the corresponding one-point cumulant in Ine. 
It must be remembered, though, that such "for- 
ward" scaling behaviour can be destroyed by the 
processes of translational averaging as well as the 
experimental "backward" box summation JlSf . 
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5. Discussion 

We have shown that features of the analytical 
solution for cuniulants in Ine can be preserved 
beyond the complication of translational invari- 
ance, and in the process elucidated the interplay 
between the same-lineage and splitting cumulants 
generated at each cascade splitting on the one 
hand, and the geometrical features on the other. 
We are, of course, tempted to apply two-point 
cumulants of Ine directly to the experimental 
energy dissipation field deduced from hot-wire 
time series and to study possible dependences 
on the Reynolds number and the flow configu- 
ration. This may be done in the spirit of naive 
discovery. We do think, however, that studies of 
different model assumptions such as continuous 
multiplicative cascade processes 1 14 , hierarchical 
shell models effects of finite inertial range 
etc. should sensibly be undertaken before taking 
the comparison with data too seriously. 
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